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ABSTRACT 

An AAT/2dF Spectrograph Survey of low Galactic latitudes targeting the putative Canis Major 
dwarf galaxy, and the (possibly) associated tidal debris of stars known as the Monoceros Ring, 
covering Galactic coordinates 231.5° <l < 247.5° and -11.8° <b<-3.8°, has revealed the 
presence of the Monoceros Ring in the background of the Canis Major dwarf galaxy. This 
detection resides at a galactocentric distance of ~18.9±0.3kpc (13.5±0.3kpc heliocentric), 
exhibiting a velocity of ~132.8±1.3 kins" 1 with a dispersion of ~22.7±1.7 kms -1 ; both of 
these comparable with previous measurements of the Monoceros Ring in nearby fields. This 
detection highlights the increasing complexity of structure being revealed in recent surveys of 
the Milky Way thick disk and Halo. 
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1 INTRODUCTION 



Galaxy formation based on a ACDM cosmology predicts larger 
structures being for med from the acc umulation of small er sys- 
tems (e.g.|Searle & Zinr]|l978tlWhite & Reeslll978tlWhitell978t 
lAbadi. Navarro. Steinmetz. & Eke 2003a b). One of the core out- 
comes of this model is that the halos of galaxies should be strewn 
with the debris from all these minor mergers. The first discov- 
ery of such a merge r within our own Milky Way, the Sagittar- 
ius dwarf galaxy llba ta. Gilmore. & Irwirll9 94). demonstrated that 
such mergers do occur and are in fact ongoing. 

Insight into the complex nature of the Galactic Halo has 
been in part accomplished by recent all sky surveys, Sloan Dig- 
ital Sky Survey (SDSS) and the Two Micron All Sky Survey 
(2MASS). Investigatin g an o verdensity of F-turnoff stars in the 
SPSS. iNewberg et alJ 1 2002) discovered the presence of a stream 
of stars suggestive of an equatorial accretion event around the 
Milky Way. Unlike the Sagittarius dwarf galaxy which is on a polar 
orbit, this accretion event may have implications on the formation 
of key Galactic structures, such as th e thick disk. The e quatorial 
orbit of the stream was investigated by ICrane et"al]l2003l) who tar- 
geted 2MASS selected M-giant stars in the range / = 150° - 240°, 
although primarily in the Northern hemisphere, finding that a sim- 
ple circular orbit model with a galactocentric radius of 1 8 kpc and 



velocity of ~220kms _1 best fit the data. The velocity dispersion 
of their sample was 20±4kms~ 1 which excludes it from being 
of Galactic origins. As more detections of this structure were made 
jYannv et alj2004Hlbata et all2003tlRocha-Pinto et all2003l) it be- 
came clear that, as proposed by Newbere et al. 1 2002), this was in- 
deed another on-going accretion event. The s earch for the progen- 
itor of the stream in the 2MASS catalog led iMartin et ail 12004a ) 
to uncover the Canis Major dwarf galaxy (l,b) = ~(240°,-8°). 

Because of the low density of stream material on the sky 
it has been necessary to employ the use of Wide-Field Cam- 
eras to look deeply into Milky Way's thick disk and Halo, not 
only confirming and enhancing previous detections but probing 
new regions leading to more detections of the Mon oceros Ring 
and the Canis Major dwarf being mad e tMartinetalJ l2004a ; 
Bellaz zini et al J 120041 : IConn et al. 112005k iMartinez-Delgad oetal] 
2005). Detectio ns of an additional structure in T riangulum- 
Andromedae iRocha-Pinto et alll2004t LMaiewski et alj|2004h are 
revealing the increasing complexity of accretion events in the outer 
regions of the Galaxy. Current debate on these detections centres on 
whether these represent separate accretion events or are the product 
of a single ongoing merger. The location of the progenitor is still 
contentious but this paper will refer to the Canis Major overdensity 
as the Canis Maj or dwarf galaxy f CMa ), following the concl usions 
of iMartinez-Delgado et aU f2005) and Mart in et alJl2005d) . 
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While the positions of the streams on the sky are highly ac- 
curate, the distance estimates to the stars are less so, requiring 
the need for radial velocity kinematic measurements to increase 
and characterize the known properties of the streams. Distances 
are th en estimated th rough photometric parallaxes as described in 
Marti n et alJ |2004b). To date, several kinematic surveys of the 
Monoceros Stream have been undertaken using the data of 2MASS 
and S PSS jYannv et alJ2004ICrane et alj2003tlRocha-Pinto et alJ 
12004 . These have confirmed a velocity gradient with Galactic lon- 
gitude and also a low velocity dispersion of o~( Vr ) ~20kms -1 for 
the Monoceros Ring. 

While most models of equatorial accretion onto disk-like 

I n 1 1 1 i t I 

Penarrubia et al. 2005; Martin et al. 2005a) predict 

the presence of multiply wrapped streams, no detections provide 
conclusive evidence supporting this hypothesis. This is most likely 
due to the pencil beam nature of most fields with regard the en- 
tire structure on the sky and given that confirmed detections of 
the stream are limited to galactic longitudes (120° < / < 240°) 
and significant areas of this structure are yet to be sampled. This 
new detection presented here is consistent with multiply wrapped 
streams. 



2 OBSERVATIONS AND REDUCTION 

In April 2004, a survey was undertaken with the 2dF (Two-degree 
Field) multi-fibre spectrograph on the Anglo-Australian Telescope 
(AAT) with two aims: firstly, to determine the kinematics of the Ca- 
nis Major dwarf, and secondly to try and locate the stars belonging 
to the Monoceros Ring through kinematic constraints. Out of ~15 
fields, eight were dedicated to this purpose, four focusing solely on 
the Canis Major region. This paper will focus primarily on results 
pertaining to the detection of the Monoceros Rin g in the C anis Ma- 
jor dwarf region, the same area studied by iMartin et alJ l2005al) ; 
hereafter Paper I. 

The 2dF instrument is capable of measuring velocities of 
~400 stars simultaneously within a two degree field of view. The 
field is divided between two spectrographs, the first using a 1200V 
grating (4600-5600A at lA/pixel); the second using a 1200R grat- 
ing (8000-9000 A at 1 A/pixel). The configuration on the second 
spectrograph was chosen to target Red Clump stars in the distance 
range 5-8 kpc, these stars fall outside the distance range expected 
for the Monoceros Ring and thus are not discussed in this paper. 
(Paper I discusses the results of the survey in the region <10 kpc.) 

The Red Giant Branch (RGB) stars in this survey were chosen 
from Sample A of Martin et al. ( 2004a). Using the photometric par- 
allax technique of iMaiewski. Skrutskie. Weinberg. & OstheimeJ 
(2003) to determine the distances, the stars in the correct distance 
range for the Monoceros Ring can then be easily selected. Deter- 
mining the exact error on the distance using this technique is dif- 
ficult and so for these observations a conservative estimate of ±1 
kpc for these stars has been assumed. The J,H,K magnitudes from 
the 2MASS catalo gue for these stars have been de -reddened using 
the dust maps of Schlege^ Finkbeine^ & Davis I fl998) and the 
asymptotic correction of Bonifacio et al. (2000). 

2.1 Reducing 2dF spectra 

Data reduction of the 2dF spectra involves using the 2dFDR pack- 
age fTavloretalll99d) which applies the flat field correction, sky 
subtraction and extracts the individual spectra from each fibre off 
the CCD. Due to an asymmetry in the Line Spread Function (LSF) 



of the spectra a custom-made reduction pipeline was constructed 
to minimise the errors this produced when comparing with radial 
velocity standar d stars; a detailed presentation of the pipeline can 
be found in Martin et al. (2005b). This pipeline limits the system- 
atic offsets in the radial velocity due to the LSF to ±5kms _1 , 
this value is not included in the errors stated. Determining the ra- 
dial velocities, v r , using several artificial standard star templates 
and taking a weighted average of the solutions avoids any system- 
atic of fset that may result f rom a peculiarity in one of the tem- 
plates iMartin et alj2005albl) . Stars then with a v > 5 kms" 1 (See 
Equation 2 of Paper I) are removed from the sample to avoid con- 
taminating the results with poor determination of radial velocity. 
Dust extinction in this region typically is E(B-V)<0.4 (taken from 
ISchlegel. Finkbeiner. & Davis I I 1998 > and only the field closest to 
the plane has higher extinction. There has been no significant im- 
pact on the 2MASS colours due to extinction. 



2.2 Determining the Fit and the Associated Error 

For each velocity and distance distribution, a Gaussian has been 
fitted to the profile using an amoeba routine IPress et alll 19921) . 
obtaining a characteristic velocity, velocity dispersion, distance and 
distance dispersion which best represents the data. The errors on the 
fitted parameters were determined by bootstrap re-sampling and 
then refitting the data. The re-sampling was undertaken using the 
errors of each quantity namely the conservative ±1 kpc for the 
distance to each star and the root-mean-square value of the disper- 
sion. The errors cited then, are the standard deviations of the fit 
determined using this method. There is no appreciable difference 
between the galactocentric and heliocentric errors in this field. The 
velocity errors also have an additional ±5kms _1 systematic error 
on top of the stated errors. 

3 RESULTS 

3.1 The Monoceros Ring behind the Canis Major dwarf 

The data in Figure Q was taken from six fields of the 2dF spectro- 
graphic survey and shows a population of ~38 RGB stars behind 
the Canis Major overdensity. These stars are represented in the solid 
histogram while the dashed histogram is the remaining stars in that 
region. The stars have been selected using the following criteria, 
231.5°<Z<247.5° and -11.8°<b<-3.8°; 12kpc<D Q <20kpc and 
a colour cut of 0.9 <(J-K) D < 1.3; 0.561(J-K) o +0.22 < (J-H) < 
0.561(J-K) o +0.36; 13.0; with v r > 0km s _1 . This overden- 
sity shows up very clearly in the galactocentric distance histogram 
and given the distance range of these stars makes them different 
from the foreground Canis Major stars as presented in Paper I. They 
form a coherent velocity and distance structure similar to that of 
previous Monoceros Ring detections. FigureQshows the 2dF Data 
with a Gaussian fit of v r = 132.8± 1.3 kms -1 and internal disper- 
sion of 22. 7± 1. 7 kms" 1 , which is comparable to the dispersion 
as measured by lYannv et alJ 12004) from SDSS data. The galac- 
tocentric distance estimates to this feature are 18.9±0.3 kpc with 
a dispersion of 0.8±0.4 kpc. See £|2.2l or Paper I for a more com- 
plete breakdown on the errors, however note again that the errors 
have an additional systematic offset of ±5 kms" 1 . To test the con- 
tamination of local dwarfs which fall inside the selection criteria 
for the sample stars due to poor extinction correction, the data has 
been tested with a more stringent cutoff of (J-K) o >0.95. This still 



A radial velocity survey of low Galactic latitude structures: II. 



3 




Figure 1. Histogram of the 2dF RGB data in the range 231.5°<1<247.5° and -1 1.8°<b<-3.8°; with a heliocentric distance cut of 12 kpc< D e <20kpc and 
the following additional cuts of 0.9 <(J-K) D < 1.3; 0.561(J-K) o + 0.22 < (J-H) D < 0.561(J-K) o + 0.36; K D ^ 13.0 and v r > Okms -1 . This field constitutes 
the region behind the Canis Major overdensity revealing a detection of the Monoceros Ring at a galactocentric distance of 18.9±0.3 kpc and a dispersion of 
0.8±0.4 kpc (right panel) with a heliocentric distance of 13.5±0.3 kpc with similar dispersion (middle panel). The velocity profile fitted to this population 
is v r = 132.8±1.3 kms -1 with an internal dispersion of 22.7±1.7 kms" 1 (left panel), the error is accompanied by an additional systematic offset of 
±5 km s _1 . The dashed lines shows the remaining data outside the distance and colour cut. For a full outline of the errors see 32.2l and Paper I. 



reproduces the detection presented here suggesting the contamina- 
tion level is low. 

Although, as can be seen in FigureQ there is no definitive ev- 
idence for a gap between the Canis Major detection and the Mono- 
ceros Ring detection this is most likely due to the large errors in the 
distance rather than the structure being somewhat continuous out 
to large galactocentric distances. The conservatively estimated er- 
ror on the distance is ~ 1 kpc, which will have the effect of blurring 
any intrinsic gap between the two structures. Contamination from 
local dwarfs although believed to be minimal will also impact on 
the ability to resolve the stream out cleanly from the Canis Major 
population. It should also be noted that depending on the formation 
scenario of the stream the stars may or may not be neatly confined 
at one distance but rather diffused over a larger volume. Since the 
origins and evolution of both structures is still being investigated, 
explaining the lack of a gap in distance is at best speculative. With 
these stars being selected from the 2MASS catalogue a measure 
of the completeness of the sample becomes possible. Applying the 
colour cuts used on the 2dF data to the 2MASS catalogue reveals 
399 stars in this region of which 364 have been observed making 
this sample ~91% complete. 

In Paper I, a northern comparison field at (240.0°, +8.8°) is 
presented for the Canis Major detection. Unfortunately, the number 
of RGB stars sampled out to large Galactic radii is too low to form 
a useful comparison with the detection presented in this paper. 

3.2 The Besancon Synthetic Galaxy Model 

Understanding what the Galactic contribution is in any given re- 
gion of sky is critical when attempting to interpret the results of 
this 2dF survey. To investigate this, th e synthetic galaxy model of 
Robin rRevle. Derriere & Picaudl2003l) is used to find out what the 
model predicts our Galactic contamination should be. This is done 
by selecting the exact same regions of sky as observed in the 2dF 
survey. Given the tight constraints in colour and magnitude used 
to select the stars in our survey (see Paper I) applying the same 
constraints to the synthetic Galaxy model to provide an idea of the 



contamination we should expect was undertaken. In particular these 
are: 0.9 < (J-K) D < 1.3; 0.561(J-K) o +0.22 < (J-H) Q < 0.561(J-K) o 
+ 0.36; K < 13.0; v r > 0km s" 1 . The cuts were applied to the 
model which extended to a heliocentric distance of 100 kpc . 

To remove the presence of noise in the model, each field was 
investigated using the online model selecting the small field setting 
with a solid angle of 1000 deg 2 . This has the effect of providing a 
smoother distribution of the stars with the properties determined by 
the centre of the field, rather than the large field option which takes 
into account changes in (l,b)° . Each field is then rescaled back to 
the size of the 2dF fields and rescaled to match the completeness 
of the 2dF sample with regard the 2MASS catalogue. These stars 
are presented in Figure [2] and in general have a much broader dis- 
tribution than could be attributed to any tidal debris with velocities 
which typically peak well away from the 2dF detections. As ex- 
pected, in the distance range of the MRi there is no thin disk popu- 
lation and with the thick disk only having ~ 15 stars in this location; 
halo stars are expected to only have a very small presence in this 
field and are not included in these results. The comparisons with 
the Besancon model provide good support for the validity of this 
detection with the usual precautions. 

The approximately 15 stars which fall in the region of the 
MRi in the model have the following properties, a velocity profile 
of v r = 143.8 kms -1 and dispersion of 28.2 kms -1 ; galactocen- 
tric distance profile of 19.7 kpc with a dispersion of 1.9 kpc. A 
Kolmogorov-Smirnoff test of the MRi stars from the 2dF sample 
and the Besancon stars finds a 4% chance of the two being drawn 
from the same distribution. This is a relatively high probability of 
coincidence however a few extra factors need to be taken into ac- 
count. While the model contains thick disk stars with an additional 
warp component, the data contains thick disk stars, CMa stars and 
MRi stars with a warp component mixed in. What effect does this 
have on our comparison? Since a completeness factor has been de- 
rived from the 2MASS catalogue and subsequently used to scale 
the Besancon model and given the 2MASS stars contains a signif- 
icant proportion of CMa stars then the numbers of thick disk stars 
must be being overestimated. This has the immediate effect of low- 
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Figure 2. Histogram of Besancon data covering the same region of sky as the 2dF dataset; no distance cuts have been applied less than the 100 kpc cutoff of 
the model. The distances are galactocentric to allow easy comparison with the 2dF data as shown in Figure^] The top panels show the thin disk population 
(stars in the age range 5-7 Gyrs) in this field fitted with a Gaussian profile at a distance of ~ 11.1 kpc and a dispersion of ~ 1.7 kpc. The velocity profile of the 
thin disk stars is v r = ~71.4 kms -1 with a dispersion of ~32.3 kms -1 . In the bottom panels are the thick disk stars (noting that these stars have a different 
density profile than the thin disk stars) they too can be fitted with a Gaussian profile and reside at a galactocentric distance of 13.1 kpc with a dispersion 
of 5.0 kpc, the velocity and dispersion are found to be v r = 1 15.9 kms -1 and 47.7 kms - 1 respectively. The smaller Gaussians represent those 15 stars in 
the same distance range as the MRi stars in FigurefTlwith v r = 143.8kms -1 and a dispersion of 28.2 kms -1 ; galactocentric distance of 19.7kpc with a 
dispersion of 1.9 kpc. 



ering the number of stars that can be expected in the distance range 
of the MRi. The true number of Galactic thick disk stars would be 
much less than the ~15 stars presented in Figure[2] especially since 
most of the stars as presented in Paper I seem to belong to the CMa 
structure and that confusion over whether they represent the warp 
seems to have been sufficently discounted. 



4 DISCUSSION & CONCLUSION 

Comparisons with the Besancon synthetic galaxy model (Figure|2j 
reveal that the peak velocity and dispersion observed in this field 
cannot be associated with any known Galactic component. The 
thick disk component in the same region as the MRi has a higher 
velocity and a broader velocity dispersion which discounts it as the 
cause of the overdensity measured in this field (see Figure^. This 
is supported by the finding that the numbers of thick disk stars in 
the model are being overestimated (after scaling) due to the pres- 
ence of the CMa structure boosting the numbers of stars in the 
region. Although difficult to quantify it is expected that the num- 
bers of Galactic thick disk stars should be much less than the ~15 
stars as currently predicted. Investigating any similarities with pre- 
vious detections of the Monoceros Ring finds encoura ging agree- 
men t in the velocity d ispersion estimates of lYannv et al.l 1 2004 ) 
and ICrane et"ai] 1200 3l) . which have values of 13 - 24 km s -1 and 
20±4 kms -1 , respectively compared to 22.7±1.7 kms -1 from 
this paper. 

The be st model of the MRi with which to compare our detec- 
tion is the Crane et al. 1 2003) circular orbit model since we target 
similar stars (M giant stars extracted from 2MASS) and since they 
reach a Galactic longitude of 1~220 in the Southern hemisphere 
and 1~240 degrees in the Northern hemisphere. Moreover, as has 
already been noticed, the velocity dispersion of our sample (a = 



22.7±1.7 kms 1 ) is statistically similar to the ICrane et alJ t 2003) 
value (20±4 kms -1 ). It can be seen on Figure 13] where our tar- 
gets are plotted as filled circles, that the velocity of our sample 
is only slightly higher than that expected by the circular model. 
Such a small discrepancy is not unexpected since the MRi pop- 
ulation certainly does not follow a perfectly circular orbit (as is 
hinted by detections at different galactocentric distances within the 
14<Dgc<20 kpc distance range). Moreover, a closer look at Fig- 
ure 2 of ICrane et alj (2003) reveals the stars at high longitude in 
their sample (1>~220 deg) also tend to have a radial velocity 20- 
30 km s - 1 higher than the circular orbit model. Thus, our detection 
of the Monoceros Ring is consistent with the current knowledge of 
this outer disk structure. 

Equatorial accretion models predict the presence of wrapped 
tidal arms and this detection provides a tentative confirmation of 
such a scenario. A plausible alternative is that both the foreground 
Canis Major detection and the background Monoceros Ring detec- 
tion represent separate accretion events. However, such a situation 
complicates the picture with the requirement of two progenitor sys- 
tems hidden somewhere in the Milky Way. Hence, we conclude that 
the probable scenario is that the Canis Major dwarf and the Mono- 
ceros Ring represent a single accretion event and this new detection 
is of a wrapped tidal arm. 
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Figure 3. Distribution of the target M giant stars compared to previous 
detections of the MRi. Stars from our sample are shown as filled circles 
whereas SDSS detections are plotted as stars. The dashed line corresponds 
to the best circular orbit model from Crane et al. 12003) for a popula- 
tion orbiting the Milky Way at D<3c=18 kpc with a rotational velocity 
of v ro t=220 kms 
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